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We experimentally study the airflow in a collapsing cavity created by the impact of a circular disk 
on a water surface. We measure the air velocity in the collapsing neck in two ways: Directly, by 
means of employing particle image velocimetry of smoke injected into the cavity and indirectly, by 
determining the time rate of change of the volume of the cavity at pinch-off and deducing the air 
flow in the neck under the assumption that the air is incompressible. We compare our experiments 
to boundary integral simulations and show that close to the moment of pinch-off, compressibility of 
the air starts to play a crucial role in the behavior of the cavity. Finally, we measure how the air 
flow rate at pinch-off depends on the Froude number and explain the observed dependence using a 
theoretical model of the cavity collapse. 

I. INTRODUCTION 

The impact of a solid body on a water surface triggers a series of spectacular events: After a splash, if the 
impact speed is high enough, a surface cavity is formed which pinches off such that a bubble is entraineoP^. 
Right after pinch-off two strong thin jets are formed^, one shooting upwards and one shooting downwards. 

An aspect in the impact on liquids that has drawn very little attention is the influence of the accompanying 
gas phase. When we take into account the inner gas in the detaching air bubble, we find a singularity in 
the velocity of the inner gas. Assuming any finite flow rate for the gas, the velocity of the gas will diverge 
because the area that the gas has to flow through goes to zero. Nature has found a way to avoid a true 
singularity by letting compressibility limit the speed of the air, but nonetheless the air plays an important 
role in the final shape of the cavity just before pinch-oflP^, and can even reach supersonic speediP. 

The main objective of this paper is to understand what determines the gas flow rate in the case of an 
impacting disc and to obtain insight into the role of compressibility effects in the air. To this end we apply 
two different approaches: First, we perform volume measurements to determine the flow rate based on 
continuity, and second we measure the air flow directly by seeding the air with smoke under laser sheet 
illumination. We compare and extend our experiments with numerical simulations, where we use one- and 
two-phase boundary integral simulations, sometimes coupled to compressible Euler equations^, to determine 
the air flow, with and without taking the dynamics of the gas phase into account. This paper builds partially 
on a previous publication by Gekle et al^ where it was shown that the gas velocity reaches supersonic speeds 
and compressibility influences the shape of the cavity in the neck region. Our new results show how the 
flow rate depends on the disc size and impact speed, for which we derive a scaling law. We experimentally 
verify the position of the stagnation point that was found numerically before^ and we determine at which 
moment compressibility starts to play a role in the dynamics of the cavity collapse. 

We have structured this paper as follows: We first give a brief description of the experimental setup 
in Section [Tl] Section |III| explains the method of volume measurements, and the results are combined with 
numerical simulations. More specifically, we measure how the air flow rate at pinch-off depends on the Froude 
number and explain the observed dependence using a theoretical model of the cavity collapse. In Section [TV] 
we perform a direct determination of the air flow velocity by seeding the air with smoke and illuminating 
with a laser sheet. Subsequently, we compare the results with the velocities that we determined using volume 
measurements. Finally, in Section[V]we discuss in detail when and how compressibility becomes important. 

II. EXPERIMENTAL SETUP 

The experimental setup consists of a water tank with a bottom area of 50 cm by 50 cm and 100 cm in 
height. A linear motor that is located below the tank pulls a disc through the water surface at a constant 
speed. This disc is connected to the linear motor by a thin rod. The events are recorded with a Photron 



FIG. 1. Two snapshots of an experiment in which a disk with a radius of 2 cm hits the water surface and moves 
down at a constant speed of 1 m/s. A surface cavity is created that subsequently collapses under the influence of the 
hydrostatic pressure. Eventually, the cavity pinches off at the depth indicated by the dashed line, and a large air 
bubble is entrained. The red arrows indicate the direction of the air flow: On the left, volume is created, resulting 
in a downward air flow. On the right the bubble volume below the pinch-off depth is decreasing, and therefore air is 
pushed upwards. 



SA1.1 high speed camera at frame rates up to 20 kHz. Our main control parameter is the Froude number, 
which is defined as the square of the impact speed Uq, nondimensionalized by the disc radius i?o and the 
gravitational acceleration g: 

JJ 2 

Fr=-^- (1) 

Two snapshots of the experiment are shown in Fig.[T] The left image shows the situation right after impact, 
when the cavity is formed. A downward flow of air is required to fill in the space that is created by the 
downward moving disk and the expanding cavity. On the right a later stage in time is shown, some moments 
before the pinch-off. Here, there is a competition between the downward moving disc and the expanding 
part of the cavity on the one hand, and the collapsing part, i.e., the region above the maximum, on the 
other. The former tends to increase the cavity volume below the pinch-off depth (dashed line), whereas 
the latter decreases it. We always observe that close to pinch-off, the violent collapse is dominant and the 
bubble volume below the pinch-off point decreases, pushing air out through the neck. As the neck becomes 
thinner towards the moment of pinch-off, the gas speed increases rapidly. The remaining part of this chapter 
is devoted to measuring this air flow and comparing the results with numerical simulations. 



III. GEOMETRIC APPROACH 



The first way in which we will quantify the air flow through the neck of the cavity is an indirect one: We 
will measure the time evolution of the volume of the cavity below the pinch-off point and calculate its first 
derivative with respect to time. This will be identified with the air flow rate through the neck. This involves 
the following assumptions: (i) The air flow is incompressible, (ii) the air flow profile is one-dimensional (i.e., 
a plug flow) and only directed in the vertical direction, and (iii) the cavity shape is axisymmetric. The first 
assumption is only violated close to the moment of pinch-off, when the air speed diverges. Compressibility 
effects at this stage are investigated iiP and its effects will be discussed in section [V] We will justify the 
second assumption partially by visualizing the air flow inside the cavity and measuring the velocity directly; 
in addition it is known from two-fluid boundary-integral simulations that the flow profile is very close to 
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FIG. 2. (a) The volume of the cavity below the pinch-off depth (dashed line) is determined by tracing the boundary 
(red line) and assuming symmetry around the central axis, (b) The volume decreases as the neck becomes thinner 
until the cavity closes, (c) After pinch-off a downward jet enters into the entrapped bubble, and the bubble shows 
volume-oscillations and cavity ripples. 
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FIG. 3. Volume below the pinch-off depth as a function of time (blue dots), determined from an experiment with 
Fr = 5.f . The vertical dashed line indicates the moment of pinch-off. Close to pinch-off the volume decrease is well 
approximated by a linear fit (green line), after pinch-off the bubble oscillates with its resonance frequency (red line: 
fit with sine function). The steady growth in volume after the pinch-off is caused by the jet entering the bubble, 
which in our data analysis is not subtracted from the measured bubble volume, see main text. 

one dimensionaP. The third assumption only breaks down in the neck-region very close to pinch-off because 
very small disturbances are remembered during the collapse I2HS1 Here, this effect is only relevant locally 
on a very small scale and can therefore be neglected on the large scale where we measure the volume. 

A. Cavity volume 

We measure the volume of the cavity below the pinch-off depth as illustrated in Fig. [2j By tracing the 
contour for every frame in a movie and invoking axisymmetry we are able to determine the volume of the 
bubble below the pinch-off depth as a function of time. Such a measurement for a disk with radius 20 mm 
and impact speed of I m/s is shown in Fig. [3j In the beginning (t < 0.022 s) the volume is increasing 
(positive slope), which means that the air at the pinch-off depth is flowing downwards. At the maximum 
(t sa 0.022 s) the flux through the pinch-off depth is zero, indicating a local stagnation of the flow at this 
depth. We will study this stagnation point later, in section |Vj After this maximum the volume starts to 
decrease and the flow is directed upwards. This continues until the moment of pinch-off which is indicated 
by the vertical dashed line in Fig. [3] A linear fit (green line) reveals that the flow rate is approximately 
constant towards the pinch-off moment. More precisely, the linear fit is the time rate of change of the cavity 
volume at pinch-off, which is equal to minus the maximum value of volume-based flow rate, <£>y = —dV/dt, 
under the assumption of incompressibility of the air. We will use this maximum flow rate <I>y to compare 
the flow rates through the neck at different Froude numbers. 
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FIG. 4. Flow rate calculated from the volume changes as a function of the Froude number in a double logarithmic 
plot. Both the experimental data (black dots) and the numerical data (red diamonds) correspond to the maximum 
value of The range of experimental data is limited to Fr « 12 by the appearance of a surface seal. The blue 
line represents the fit $y = 1.23-FV 1/2 + 1.01. 

After the pinch-off there is a clear oscillation of the volume together with a slow apparent growth of the 
bubble. The growth is caused by the liquid jet that is entering the bubble (Fig. [2j;), as the amount of air 
is fixed after the pinch-off. Since our focus is on the behavior before pinch-off, we chose not to correct 
the bubble volume in the image analysis by subtracting this jet volume. Also, making such a correction 
would be complicated by the fact that the jet is imaged through the refracting, curved interface of the air 
bubble. Nevertheless, we determined the frequency of the oscillation by fitting a sine function (red line) 
after correcting for the slightly positive slope. For the conditions of Fig.[3]the measured frequency is 143 Hz, 
which is close to the expected resonance -or Minnaert- frequency^ of 138 Hz^. This agreement was also 
noted by22D for the impact of freely falling objects in water. 



B. Air flow rate 

A characteristic quantity concerning the gas dynamics in a collapsing cavity is the air flow rate, defined 
as the volume of air that is being displaced per unit time close to pinch-off. From Fig. [3] we infer that in 
approach of the pinch-off point this flow rate becomes constant and can be determined as the maximum 
slope of the volume as a function of time (green line, Fig. [3]), i.e., the air flow rate through the neck equals 
the rate of change of the volume of the cavity below the pinch-off depth, of course under the assumption 
that the gas flow remains in the incompressible limit. We subsequently non-dimensionalize this air flow rate 
$y = —dV/dt with the disk radius and the impact speed, namely, <&y = ^v/(RqUo), where the asterisk 
denotes a dimensionless value. We determined the flow rate for a number of different disk radii (ranging from 
15 to 30 mm) and impact speeds (0.45-1.30 m/s), the results of which are shown in Fig. |4]where we plot the 
dimensionless flow rate <&y versus the Froude number Fr on a double-logarithmic scale (black dots). The 
experimental range is limited by the appearance of a surface seal at high impact speeds, where the crown 
splash is pulled inwards due to the air flow induced by the disc and closes the cavity at the surface. This 
surface seal usually has a significant influence on the cavity shape and dynamics^! as well as the gas flow 
rate in the neck, so all of the experiments reported here are without surface seal. Although the experimental 
data foll ows an apparent power-law, extending the experimental range in Fr with our boundary integral 
cod d 8 * 19 ! reveals that the results do not lie on a straight line (Fig. [4j open red diamonds) . This suggests that 
there does not exist a pure power-law. 

An analytical argument - Using the assumption that the cavity expansion and collapse take place in 
horizontal non-interacting layers of fluid -an assumption that was successfully used in Bergmann et aZP^L 
we will now shed light onto the behavior of the air flow rate through the neck as a function of the Froude 
number. We will provide an analytical argument in this subsection, which we work out in detail based on 
the model of Bergmann et aZP^in Appendix [A] For convenience from hereon we will take the z to mean the 



5 




FIG. 5. A schematic view of the setup. A laser sheet shines from above on the disk, illuminating the interior of the 
cavity after the disk has impacted the water surface. We insert smoke in the top part of the container and when the 
linear motor pulls the disk through the water surface at a constant speed, the smoke is entrained into the cavity. 

depth below the undisturbed water surface, i.e., z = at the latter and increases with depth. 

The quantity that we aim to calculate is the time rate of change of the cavity volume V = dV/dt, i.e., 

dV d [ ZA ™ , , . l2 , 
* vS = n[r(z,t)] 2 dz, (2) 

where it is understood that the expression needs to be evaluated at the pinch-off time. Here, r(z,t) is the 
cavity profile, z^ sc {t) is the vertical position of the disc and z c is the pinch-off depth. We obtain 

$> v = - I 2irr(z,t)r(z,t,)dz - -kRqUq, (3) 

J z c 

where the last term is due to the downward moving disc and r = dr/dt denotes the radial velocity of the 
cavity wall. 

To approximate the integral in Eq. ^ we subdivide the expanding and collapsing cavity -at times close 
to the collapse- into regions that we can express analytically using the model provided by Bergmann et 
al™. This is done in Appendix [a] and leads to 



or $y = A Fr 1 / 2 + B , (4) 

with A and B numerical constants. To test this relation we extended the experiments of Fig.Hby performing 
boundary integral numerical simulations in order to cover a wide range of Froude numbers^T The obtained 
results are added to Fig. [4] using red diamonds. There is a good agreement with the experimental data, and 
the non-constant slope is clearly visible. A fit to the simulation data confirms Eq. Q and gives A sa 1.23 
and B 1.01. 
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IV. FLOW VISUALIZATION 



In this Section we perform a direct determination of the air flow velocity by seeding the air with smoke 
and illuminating with a laser sheet (see also Gekle et altty , the results of which we will subsequently compare 
to the velocities that were determined indirectly and independently using volume measurements. We will 
first describe the method and results of the flow visualization that we used to measure the air flow inside the 



FIG. 6. A snapshot of the cavity with an overlay of a recording of the illuminated smoke. The smoke particles are 
artificially colored orange in this figure. The size and position of the employed correlation window is indicated by 
the yellow square. 



cavity. Before doing the impact experiment we fill the atmosphere above the water surface with small smoke 
particles. When subsequently the disk is moved down through the water surface, the smoke is dragged along, 
and fills the cavity created below the surface. We illuminate a thin sheet of the smoke using a 1500 mW 
diode laser line generator (Magnum II) and record the experiment at a recording rate up to 15 kHz by placing 
the high speed camera perpendicular to the laser sheet (Fig. [5]). The smoke consists of small glycerine-based 
droplets (diameter ~ 3 /im), produced by a commercially available smoke machine built for light effects in 
discotheques. A simple analysis shows that the particles are light enough to neglect all inertial effects at 
least in the range of accelerations that we can measure experimentally: At a velocity difference of 10 m/s 
the Reynolds number is ~ 2, meaning that we can assume Stokes drag. Knowing the force on the particle 
as a function of the velocity difference and the mass of the particles, we can calculate the movement of the 
droplets in an accelerating flow. We find that the particles follow the flow up to 25 m/s with a velocity lag 
less than 2%. 



Correlation technique - We determine the speed of the air in the neck by applying an image correlation 
velocimetry (ICV) technique^. ICV differs from Particle Image Velocimetry (PIV) in the sense that we 
do not resolve discrete particles in our images, but we correlate smoke patterns instead of smoke particles. 
Figure [6] shows the cavity with the illuminated smoke as an overlay, where the smoke is colored orange 
artificially for clarity. The actual measurements are done on a closer view of the cavity. The correlation 
is performed on a square correlation window, indicated by the yellow square. The width of the correlation 
window is 160 pixels, corresponding to 8.8 mm. In the latest stages we switch to a correlation window of 96 
pixels (5.3 mm) wide, anticipating for the smaller neck radius. The measurements are insensitive to small 
changes in the shape, size or position of the correlation window. The size of the window is optimized for 
quality of the cross correlation. 

To improve the ICV analysis we used an image subtraction technique that we describe in Appendix [B] A 
subtraction technique similar to the one that we use here has been used previously for double-frame PIV 
images^!, where it was found that if the displacement of the particles is too small between a pair of images, 
the displacement peak in the correlation is biased. This bias is related to the particle size in pixels and the 
displacement in pixels. In our case this length scale does not exist because we do not resolve separate smoke 
particles in our experimental setup. Instead of calculating the expected bias, we identify biased values by 
their departure from the global trend of the data (Fig. 7j inset). As a remedy for the bias, we artificially 
increase the displacement by skipping frames. The smaller the velocity, the larger the number of frames we 
skip. In addition to this we note that the bias is less pronounced compared to the case irP^l because we 
construct the image pair from three images instead of two. 

The biased data and other spurious data is removed by making an objective selection based on the peak- 
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FIG. 7. The vertical air velocity through the neck as a function of the neck radius R, measured in an experiment 
with Fr = 5.1 in three different ways: (i) Directly, using smoke particles (diamonds), (ii) indirectly, using a smoothing 
polynomial fit to bubble volume of Fig. [3] (blue line) , and (iii) indirectly, using a constant flow rate approximation, 
determined at pinch-off (cf . Fig [3j black line) . The different colors of the diamonds correspond to different numbers 
of frames that are skipped in the cross-correlation (see main text). The inset shows the same vertical velocity data 
measured using method (i) for two different values for the peak-to-peak ratio A: For A > 1.5 (orange dots) we find 
strongly biased data, which are eliminated using a higher threshold (A > 3.5, black dots). 



to-peak ratio of the correlation. This ratio is defined as the ratio between the two highest peaks in the 
correlation: A = p\/p2- The inset of Figure [7] shows the effectiveness of this selection method. We set A 
to values between 3.5 and 5.0, depending on the specific measurement, so that almost all spurious data is 
removed. Taking higher values for A removes too many valid data points; lower values allow for too many 
biased data points. 

In Fig. [7] we compare the air speed that we measured directly using smoke particles with the air velocity 
that we calculated indirectly using the change in volume of the cavity, as discussed in the previous Section. 
The air speeds are plotted versus the neck radius R(t) at pinch-off depth instead of time; time increases 
from right to left in the figure, i.e., towards smaller values of R. The blue line is obtained using a polynomial 
(smoothing) fit to the volume-time data of Fig. [31 determining the flow rate $y(i) from the time derivative 
of this fit [Eq. ||2j)], and finally dividing by wR(t) 2 to obtain the velocity. We find a very good agreement 
between the direct (smoke) and the indirect (volume) measurements. 

Finally, the black line in Fig. [7] is obtained by setting the flow rate to a constant value, namely to that 
corresponding to the time derivative of the volume curve just before pinch off (the green line in Fig. [3]) . We 
observe that at early times (large R) there are large deviations from the other two datasets. This stands to 
reason, since at these times we are still far away from the pinch-off moment, and the gas flow rate in the 
neck has not yet become (approximately) constant. Close to pinch off however, for R/Ro < 0.4, we find 
that the constant flow rate approximation and the smoothing fit both provide the same air speed. 



V. THE ROLE OF COMPRESSIBILITY 

The fact that the air flow rate becomes constant together with the surface area of the neck becoming van- 
ishingly small suggests that the velocity in the neck diverges towards pinch off. However, as was mentioned 
in the introduction, a real singularity of the air flow velocity is prevented by compressibility effects. In a 
previous publication we presented a directly visible effect of the compressed gas flow, namely the upwards 
motion of the position of the minimum neck radius^. This upwards motion was seen both in experiments 
and in simulations that take the compressibility into account, and is absent in simulations that neglect com- 
pressibility. In the same paper we reported that, next to this upward motion of the neck, the extremely fast 
airflow affects the smoothness of the neck. Especially this last effect is important, since it is i n cont radiction 
with the assumptions in theoretical pinch-off models where the neck is assumed to be slendejP^J. 

The question that we intend to answer in the present Section is how the effects of compressibility show 
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FIG. 8. The location of the stagnation point z a tag with respect to that of the pinch-off point z c as a function of the 
neck radius R. Note that when that the stagnation point lies below the pinch-off point, z c — z s tag is negative. Time 
increases from right to left (decreasing R) . The dots are experimental data, obtained by volume measurements of four 
different experiments, where each color corresponds to a different experiment. All experiments were performed with 
disk radius Rq = 2.0 cm and impact speed Uo = 1.0 m/s, i.e., Fr = 5.1. The green line is the result of a two-phase 
boundary integral simulation without taking compressibility into account [type (ii)]. The red line is obtained by a 
two-phase boundary integral simulation which includes a compressible gas phase [type (hi)]. Oscillations in the red 
line are a numerical artifact due to wave reflections in the compressible domain, see Gekle et alW for details. 



up in the measurement of the cavity volume and the air flow rate that can be deduced from it, as was 
presented in Section |III| of this work. More specifically we will investigate the position of the stagnation 
point of the flow in the cavity (see below) and the air flow rate towards the pinch-off moment. Following 
the method we used irP, we will compare our experimental results with three different types of boundary 
integral simulations: (i) a single phase version, in which only the water phase is resolved, (ii) a two-phase 
version where both the liquid and the gas flow are resolved as incompressible inviscid media, and (iii) a 
compressible gas version where the compressibility of the gas phase is taken into account by substituting the 
incompressible axisymmetric gas phase equations by one-dimensional compressible Euler equations at that 
moment during the collapse when compressibility effects start to become significant. More details about the 
numerical method can be found irpm. 



Stagnation point - Just above the disc the air must move downwards at approximately the same speed as 
the disc, whereas simultaneously, towards closure, the air in the neck is moving upwards. This implies that 
somewhere in between there will be a stagnation point. We will estimate the location of this stagnation 
point as follows: The first step is to extend the analysis of Section III where we tracked the volume below 
the pinch-off depth in time, to any depth z below the pinch-off point. For every depth z this will provide 
us with a curve similar to that in Fig. [3] and by determining the time coordinate of the maximum we find 
the time i s tag at which the averaged^ flow rate (~ V) at that depth z = z sta , g is zero. This point we then 
interpret as the location of the stagnation point z s tag(istag), which involves the assumption that close to the 
pinch-off moment the flow in the neck region becomes predominantly homogeneous and vertical. In Fig. [8] 
we plot the measured location of £ s t ag for three different realizations of an experiment with a radius of 2 cm 
and an impact speed of 1 m/s. When we compare the experiments to a two-phase incompressible boundary 
integral simulation [type (ii)] (green line in Fig. [8]), we find a considerable discrepancy between the two for 
small values of the neck radius R. If we however use the compressible version of the simulation [type (iii)], 
the agreement becomes much better (red line in Fig. [8| see also Gekle et a/.^) , confirming the importance 
of compressibility in this limit. Note that experiments and both simulations do converge for larger values 
of R, where compressibility effects play no role. 

The agreement is not perfect, however, which can partly be traced back to the technical difficulty of 
obtaining reliable values for z stag from the experiment (which reflects in the large spread between the three 
different realizations) and partly to the fact that its determination neglects compressibility in a subtle way: 
Although in the experimental data compressibility is of course necessarily reflected in the shape of the 
cavity, the method of obtaining the air flow rate from it (namely by determining the time rate of change 
of the cavity volume) neglects compressibility in the air phase. A difference with the actual location of the 
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FIG. 9. (a) The dimensionless derived air flow rate &y = $v / (RoUo) (from the time rate of change of the cavity 
volume) as a function of the dimensionless neck radius R/Ro in an impact experiment with disc radius Ro = 2 cm 
and impact speed Uo = 1 m/s (Fr = 5.1). The black dots represent experimental data. The red line is obtained 
using a one-phase simulation [type(i)], which excludes the air phase. The green line is a two-phase boundary integral 
simulation without compressibility [type (ii)]. Finally, the blue line is the result of a two-phase boundary integral 
simulation which includes a compressible gas phase [type (iii)]. (b) Comparison of the dimensionless derived air flow 
rate & v [blue line; the same curve as in (a)] and the true air flow rate both plotted versus R/Ro. The two curves 
diverge from each other below R/Ro ~ 0.2. 



stagnation point is therefore expected for high gas velocities (i.e., small neck radii). 

Air flow rate - As explained in the introduction of this Section, we can compare the air flow rate in 



the neck in experiment and simulation directly, by comparing the experiment al v elocities (cf. Section IV I 
as was done irP, but also indirectly, by using the volume analysis of Section III both in experiment and 



simulation. This second method, the results of which will be presented now, enables us to distinguish the 
effect the compressibility of the air has on the cavity wall (which is included in the analysis) from the pure 
compressibility of the flow (which is not included) . 

To do so, it is convenient to from now on distinguish the true air flow rate from the derived air flow 
rate, i.e., the one obtained from the time rate of change of the cavity volume. In Fig. [9^a) we plot the 
non-dimensionalized experimental derived air flow rate in the neck, 

v ~ R 2 U ~ R*U dt ' [ j 

as a function of the dimensionless neck radius R(t)/Ro (black dots), again for Fr = 5.1. Repeating the 
experiment results in an uncertainty in the magnitude of <fry (corresponding to the spread of the experimental 
data in Fig. El), but the behavior as a function of time is always the same: The derived air flow rate in the 
neck reaches a maximum, and approaches a finite value towards the moment of pinch-off. We compare this 
result with those of the three different types of boundary integral simulations: 

The one-phase code [type (i)] predicts a steadily increasing derived air flow rate, which seems to level off 
to a constant value towards pinch off (R/Ro — > 0). This is the red line in Fig. [9|a). 

The two-phase incompressible version [type (ii)] predicts a maximum at a location which is reasonably 
comparable to the experimental one, but after that decreases toward zero at the pinch-off moment (the 
green curve in Fig. |9ja)). Since both phases are incompressible, this stands to reason: The pressure in the 
cavity rises instantly because of the divergence of the air velocity u g in the shrinking neck. This pressure 
decelerates the cavity wall, which in turn decreases the derived air flow rate, which should go to zero in 
the R/Ro — > 0-limit: In the context of incompressible flow, a finite derived air flow rate would result in an 
infinite air velocity in the neck and consequently an infinite pressure within the cavity. Here it is good to 
note that for this two-phase incompressible code the derived and true air flow rates are actually identical, 
due to the incompressibility of the air phase. 

The two-phase compressible simulation [type (iii)] also predicts a maximum for at a location similar 
to the two-phase incompressible code and the experiment, but then decreases to a finite value for R/Ro —> 0, 
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just like the experiment. Clearly, and in contrast with the other two versions of the simulation which behave 
poorly, the agreement with the experiments is qualitatively very good and quantitatively satisfactory. All 
three types of simulations and the experiments all converge for larger R/Rq ~ C(l), which is expected since 
airflow effects (let alone compressibility of the air phase) are small or even negligible in that regime. 

The final question that we want to address is the difference between the true air flow rate (which incor- 
porates all compressibility effects) and the derived one (which only includes the effects of compressibility 
on the cavity wall). In experiment it is impossible to obtain the first quantity at the required precision, 
because its determination includes measurement errors in both air velocity u g and neck radius R. The 
two-phase compressible simulation technique however does offer a way to look at this difference: In Fig. [9^6) 
we compare the derived air flow rate <&y (the same curve as the blue one in Fig. ^a)) to the true air flow 
rate 4>* (green curve), which is calculated from u g (t) and R(t) as 

= WfT = WfT™^) 11 ^ 2 > ( 6 ) 

both as a function of the dimensionless cavity radius R(t)/Ro- Clearly the two curves coincide above 
R/Ro ~ 0.2, but start to depart from one another below this value, indicating that here the compressibility 
of the air itself becomes significant, in good agreement with what we concluded from the previous plot 
(Fig. [i^a)). We observe that the true air flow rate goes to zero for R/Rq —> (and incidentally not quite 
unlike the two-phase incompressible curve (green) in Fig.^a)). This is of course what should happen, since 
the gas velocity in the neck needs to remain finite at all times. The difference between the two curves is 
related to the rate at which the gas in the cavity is compressed. 



VI. CONCLUSIONS 



We have measured the air flow inside the neck of a collapsing cavity that was created by the impact 
of a circular disc on a water surface. More specifically we have performed and compared two types of 
experiments: First, we did indirect measurements, using the time rate of change of the cavity volume as 
a measure for the air flow rate in the neck, thereby neglecting compressibility of the air inside the cavity. 
Secondly, we performed direct measurements of the velocity in the neck of the cavity using image correlation 
velocimetry. Numerical boundary integral simulations of three different types have been used to evaluate 
and discuss our experimental findings. 

For the complete experimentally available range of Froude numbers we showed that there is a very good 
agreement between the indirectly measured air flow rate and the boundary integral simulations. With the 
simulations we were able to extend the range of experimentally attainable Froude numbers, which revealed 
that the air flow rate is not a pure power-law of the Froude number. We formulated an analytical argument 
revealing that the dimensionless air flow rate should scale as AFr 1 ^ 2 + B. Such a scaling compares well with 
experiments and simulations for A « 1.23 and B m 1.01. 

By performing careful image correlation velocimetry experiments with a smoke-filled cavity we have been 
able to directly measure the air flow for relatively low air speeds, corresponding to R/Rq > 0.3. In this region 
we found excellent agreement with the gas velocities that we calculated from the indirect measurements of 
the air flow rate and the neck radius R(t). 

Due to the very high air speed close to the moment of pinch-off (R/Rq < 0.2) compressibility of the air can 
not be neglected anymore. We have demonstrated this by comparing experimental results to three types 
of numerical simulations: (i) one-phase boundary integral simulations, (ii) two-phase boundary integral 
simulations with an incompressible gas-phase, and (iii) a compressible gas version of the second type of 
simulations that include the gas phase as a compressible fluid. We analyze the time evolution of both the 
location of the stagnation point in the gas flow and the derived air flow rate and explain our experimental 
observations in terms of the three types of simulations. The main conclusion is that the behavior that 
we observe in the experiments can only be reproduced by the simulations if compressibility is taken into 
account. 



Appendix A: Derivation of the scaling law for >l>[ 



In this Appendix we derive Eq. Q, the main result of §IIIB starting from the description of the cavity 



proposed irPD. The starting point is the two-dimensional Rayleigh equation for the cavity wall r{z, t), which 
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cross 



FIG. 10. For the derivation in this Appendix, the cavity close to the pinch-off moment needs be divided into four 
regions: The expansion region (A), between the location of the disc Zdisc and the location of the maximum zm, 
where the cavity expands against hydrostatic pressure; the contraction region (B), between zm and the point z cro3S 
where the cavity reaches the disc radius again where the hydrostatic pressure approximation [Eq. (A2i] is matched 
to the inertial approximation [Eq. ( |A4[ )]; the collapse region (C) between z cross and z*, characterized by continuity; 
and the self-similarity region (D), between z* and the pinch off location z c , which is in addition characterized by a 
coupling between the vertical and horizontal coordinates. 



originates from integrating the Euler equations in uncoupled horizontal layers of flow from some far away 
point i?oo to the cavity wall 



log (r/R^ — irr) + |j 



gz 



(Al) 



in which r = dr/dt and g is the acceleration of gravity. 

This equation is solved in two different limits to describe the different regions in Fig. |10| The first one 
is to describe region A and B, taking for every depth z the moment tM{z) of maximum expansion as a 
reference point. With r(i/^) = Rm(z), r(tnj) = 0, we can neglect the second term in Eq. (Al) and replace 
the slowly varying logarithm in the first term by a constant 2 ^ j3 = log (Rm /Roo an d solve 



r(z,t) 2 =R M (z) 2 



gz 



{t-t M (z)f 



(A2) 



IrP it was shown that 



tu{z) 



TT i ^expaPexpa 

U gz 



(A3) 



in which the first term represents the time span needed to arrive at depth z and the second the amount of 
time to expand to the maximum radius. Here a expa /3 expa is a constant!^. 

The second approximate solution corresponds to the small R limit in the collapse regions C and D of 
Fig. 10 in which both the driving pressure gz and the inertial term f 2 /2 can be considered small when 



j logfr/i^oo)! ^> 1/2. This then leads to ^{rf) = which is readily solved to give: 

r(z,t) 2 = 2a ctra R U {t coU {z) - t) , 



(A4) 



in which a ctra is a constant and t co u(z) is the (virtual) closure time of the cavity at depth z. At any depth 
z the approximate solutions are tied together at the maximum (where a solution Eq. ( A2 ) with /3 = (3 expa 
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is matched t o a s olution with j3 = /3 ctra ) and at th e mo ment t cross (z) when r(z,t) — Rq again. Here, the 
solution Eq. ( A2 ) with ft — (3 ctra is matched to Eq. ( A4 1 . More details can be found irpSl 

The quantity we want to calculate is Eq. ^ty, which contains the time derivatives of Eqs. (A2| and (A4), 
which are: 



d_ 
dt 
d_ 
dt 



r(z,t) 2 ) 
r(z,i) 2 ) 



-2^{t-t M {z)) (regionA,B) 
-2a ctra RoU (regionC,D). 



(A5) 



which subsequently need to be evaluated at the m oment o f pinch-off t = t c , for which it was derived irP^ that 
it is independent of the impact speed: t c — C 2 y/ Ro I y^- Inserting this expression together with Eq. (A3) 
into the first Eq. ( A5 ) gives 

d (r(z,t c ) 2 )= (A6) 



dt 

-2 9 4 I C 2 




Uq gz 



Finally, we need to integrate the second Eq. (A5) and Eq. (A6) over z between Zdi SC and z c . This is a 
straightforward calculation which gives the following lengthy result 

«disc d 



dt 



C 2 



(r(z, t c 



Ro9 (z d 



dz = 



(A7) 



ftexpa 

2&expaRoUo( z disc ~ z m) 

2 9 



c 2 

ftctra 



M 







M 



2a ct raRoU (z c 



j +2 a exp aP expa RoUoi2M 

Pctra 

z+) - 2a c t ra R Uo(z* - Zc) ■ 



We now use that all length scales Zdisc, %M, z cross , and z. M scale as RoFr 1 ^ 2 , except for the differe nce (z* — z c ), 
which due to the self-similarity in the neck radius scales as Rq. This means that the above Eq. (A7) has the 
folowing form 

' d (r{z,t c f)dz = (A8) 



dt 



-K^RogR^Fr + k 2 ~ R^Fr 



■" ^ 3/2 ~*~ K 3 R 2 U Fr 1/2 - K4 R 2 U : 



in which K1-K4 are positive numerical constants, which depend on the a's, B's and the proportionality 
constants in the scaling laws for the length scales Zdisc, %Mi z cross , 2*, and (z* — z c ) . By writing g = 
UqRq 1 Fr^ 1 in the first two terms we finally obtain: 



dt 



{r{z,t c ) 



2 )dz 



(A9) 



(-/si + k 2 + k 3 ) R 2 Q U Fr 1/2 - k 4 R 2 U . 
If we now insert the above result in Eq. (|3| we obtain 



4n 



7r(Ki 



2ir r(z,t) r(z,t,)dz — ttRqUq 
- k 2 - k 3 ) R 2 U Fr 1/2 + tt(k 4 - 1)F%U 



(A10) 



which then leads to 



V 



AFr 1/2 + B. 



(All) 



with A = 7t(ki — k 2 — K3) and B = tt(k4 — 1) 
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Appendix B: Image subtraction technique 



Mainly due to reflections from and refraction at the free surface, there are structures visible in the correla- 
tion window that move slowly compared to the typical gas velocities that we want to measure. A correlation 
between two unprocessed images gives a strong correlation peak close to zero because these structures are 
dominating the image, and thereupon also the cross-correlation. Standard background subtraction is not 
able to remove these features since, because of their refractive and reflective nature, they appear and dis- 
appear at unpredictable instances in time and are not stationary. Instead we use the difference between 
subsequent images in the following way: We start with three images /„, I n +i, and I n +2- After applying a 
low pass filter we create two new images by mutual subtraction of the original three images: J n = I n +\ — I n 
and J n +i — I n +2 — In+i- We then apply a min-max filtei^to these images J n and J n +i, followed by the 
cross correlation of J n and J n +i- On the result of the correlation we apply a multiple peak detection to find 
the highest peak pi and the second-highest peak p2 . We determine the position of the highest peak with 
sub-pixel accuracy by a gaussian fitting routine. 
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